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Centaurus A is the nearest radio-loud AGN and is detected from radio to very high energy gamma- 
rays. Its nuclear spectral energy distribution shows two peaks, one in the far-infrared band and 
another at about 150 keV. By assuming the second peak is due to the electron synchrotron emission 
and the power index for the differential spectrum of the very high energy cosmic ray proton to be 2.7 
we show that only pp interaction is responsible for the observed GeV-TeV emission from Centaurus 
A. We also found that indeed many very high energy cosmic ray protons from Centaurus A can 
arrive on Earth thus supporting the recent observation of two events by Pierre Auger Observatory. 
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I. INTRODUCTION 
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^^ \ Centaurus A (Cen A or NGC 5128) is the nearest active radio galaxy with a distance of approximately 3.8 Mpc and 

CN ■ redshift z = 0.002. Because of its proximity, Cen A is one of the best studied radio galaxies although its bolometric 

' luminosity is not large by an AGN standard. Optically, Cen A is an elliptical galaxy undergoing late stages of a 

r^ merger event with a small spiral galaxy. From Cen A, sufficiently large amount of photometric data is available to 

fS i' build a well sampled spectral energy distribution (SED). The emission from the nucleus of Cen A has been observed 

JL \ in radio, infrared, X-ray, 7-ray[iH3 and also in the GeV-TeV range lp-lia . Recently the Fermi Large Area Telescope 

ij ' (LAT) has also observed Cen A in the energy range 0.2 to 100 GeV[12]. Observations in different wavelengths show 

(~| that Cen A has FR I morphology having two radio lobes, non-blazar source with a jet inclination of about SO'^. It 

'— '. has a central supermassive black hole of mass msH ^ (0.5 — 1.2) x 10^ Mq. The nuclear SED shows two peaks, one 

^^ ' in the far-infrared band (~ 4 x \0~^eV) and another at about 150 keV. 

►^ High energy 7- Rays are produced due to non-thermal process where the particles are accelerated to ultra relativistic 

l/~\ ' energies. Detection of GeV-TeV photons from Cen A signifies that, Cen A has the potential to produce very high 
l/~j , energy cosmic rays. Recently Pierre Auger Collaboration reported a correlation between ultra high energy cosmic 
■TJ" ■ rays (UHECRs) and nearby AGN within ~ 75 Mpc. In particular two of these events fall within 3.1° around Cen 
^^ ' A, thus strengthening the possibility that this object is a strong candidate for UHECR source. By assuming the two 
[ — . events are from Cen A the expected rate of UHE neutrinos in detectors like IceCube are estimated [Ij, [iJI . By using 
f^ ' the same hypothesis the diffuse neutrino flux from Cen A are estimatedfl^. Also the flux of high energy cosmic rays 
^^ ■ as well as the accompanying expected secondary photons and neutrinos are calculated from hadronic models [16|. 
'""' ■ It is logical to argue that the astrophysical objects which are producing UHECRs also produce high energy 7-rays 

J> , due to interaction of the UHECRs with the background [16l - l20l | . We have to also keep in mind that apart from the 
hadronic interaction [pp and p^) with the background, there are leptonic processes (e.g. electron synchrotron. Inverse 
Compton (IC) Scattering, etc.) also responsible for the production of high energy 7-rays and efficiencies of both the 
jrt ' hadronic and leptonic processes depend on the background particles and/or the magnetic field. By using the currently 
available GeV-TeV gamma-ray data from Cen A, it has been stressed in the literature that only hadronic interaction 
is insufficient to explain the two events from Cen A|2l|. 

Our paper is organized as follows: In section 2, we discuss different processes which are responsible for the production 
of gamma-ray in Cen A and how it is compared with the observed photon spectrum. Section 3 is devoted to estimate 
the number of events from Pierre Auger Observatory (PAO) and its relation with the CR normalization constant. A 
summary of the observations of very high energy gamma-rays above 250 GeV by different experiments and calculation 
of expected events from these experiments are discussed in section 4. We briefly conclude our results in section 5. 

II. GAMMA-RAY EMISSION PROCESS 

The main mechanisms responsible for production of 7-rays in astrophysical objects are the: synchrotron emission 
by charge particles (protons and electrons), inverse Compton (IC) scattering of electron by background photons, 
inelastic scattering of high energy cosmic ray proton with the background protons and the photons, and also the 
photo-de-excitation of nuclei. Depending on the projectile energy and the nature of the surrounding environment of 



the astrophysical object some or many of the above processes can be effective. Here we study the case of Cen A where 
we have to analyze which are the most probable mechanisms to produce the high energy ^-ravs[22i[23{. 

In general for mildly relativistic systems like AGNs, proton synchrotron radiation is an inefficient process in com- 
parison to electron synchrotron radiation and the energy loss rate of proton is {rnp/nie)'^ ^ 10^"^ times slower than the 
electron. Also emitted photon energy from the proton is {nip/mf,)^ '^ 6 x 10^ times smaller than the photon emitted 
by the electron of same energy as proton. The above analysis is also true for a relativistic systemi24j where proton 
synchrotron process is suppressed in comparison to the leptonic one. So in this case we neglect the proton synchrotron 
radiation process as compared to the electron synchrotron one. On the contrary in Gamma-Ray Bursts, where the 
jet is ultra-relativistic, synchrotron emission by proton can be an efficient process. Also the photo-de-excision con- 
tribution to the 7-ray is very small in Cen A environment, so we do not take this into account here. In general, the 
energetic electrons which are Fermi accelerated in the jet of Cen A can have a power spectrum given by 



dEe 




Ee > E^h 



The constant A^ has the dimension oi photon / keV / cm? / s and Eq is the energy normalization constant which can vary 
depending on the range of energy one considers. The above electrons will lose energy through synchrotron emission 
and the break in the electron energy Ei,,b is given by 

and correspondingly the break in the photon energy is given by 

,,_ 277r3eme ^r^ M V^' ^^"^ ^ 

where ctt — 6.65 x 10~^^ cm^ is the Thompson scattering cross section and 13 = [(F^ — 1)/F^]^/^. Also F, ^b , Ai"*"* and 
fes are respectively the bulk Lorentz factor within the jet, fraction of energy carried by the magnetic field, observed 
time variability and ratio of shell expansion time to synchrotron emission time which is given by 

The L°'"' corresponds to the observed luminosity in the observed range of photon energies. 

With the synchrotron self-Compton model, M. Chiaberge et al. have reproduced the whole nuclear emission of Cen 
A by takin g th e bulk Lorentz factor F < 3 — 5 and these values are consistent with the mildly relativistic [25|,[26| or a 
little more [23] proper motions observed on sub-parsec scales [28|. Also in ref.|2^ using SSC model they estimate the 
F < 2.5 in parsec scale. The observation of Cen A during 1973 to 1983 in X-rays below 100 keV has shown greater 
than an order of magnitude variation on time scales of years |30l - [33 | . Later observations have also shown variability 
on similar timescales [3j, [sj but no high luminosity state was detected since 1985. 

In the present work we assume that the second peak at 150 keV in the nuclear SED region is solely due to the 
electron synchrotron emission in the nuclear region of Cen A and we use the energy normalization Eq — 100 keV [35|. 
The free parameters in the electron synchrotron model are F, £^b , At"'"' and fes- But from the above discussions we 
know that jet in Cen A is mildly relativistic and around 100 keV range the time variability is of order years. In order 
to fix the observed photon energy E°''^ = 150 keV by synchrotron emission we take the parameters Ai°^* ~ 1 yr, 

F = 4.2, ^B = 0.2 and fes =0.1. The two parameters F = 4.2 and At"'"' ~ 1 yr imply that the jet in Cen A is mildly 
relativistic and the time variability is order of year in this keV range of the photon energy. Here we further assume 
that the same time variability also holds around 150 keV. The parameter ^b = 0.2 signifies that the equipartition 
energy in the magnetic field corresponds to 20% of the total luminosity. Using the above values of the parameters, 
the magnetic field is given by, 




'3 /A,o6s\-i/ L°^* \ 1/2 



V lyr J V3.78X lO^^erg/s 



G. (5) 



Also the time variability corresponds to a distance scale of r '-^ 2F^Ai°'"' ~ lO.Spc which is in the nuclear region of 
Cen A. 



Cen A has a dusty torus within 100 pc of the black hole, with high column density and there can be infrared 
emission from it. Additionally we also assume that the first peak at around ~ 4 x 10""' eV in the nuclear SED 
has contribution from electron synchrotron emission. But for electron synchrotron emission to be operative in this 
region the time variability is much smaller compared to the one for 150 keV peak and one can consider the minimum 
variability time scale At°^^ ^ 1 day |36[. To fix the first peak we keep the bulk Lorentz factor F — 4.2 as before 
but vary ^b and fes- In this case both ^b and fes are small compared to their corresponding values for 150 keV . 
This shows that for low energy peak the synchrotron cooling is very very slow and also the equipartition energy in 
the magnetic field is very small. As discussed in ref. [24| , the energy dependence of the differential photon spectrum 
for electron synchrotron emission is proportional to E-y below the break energy ^--^ 4 x 10"** eV. In a very recent 
paper by Araya and Cui'37] it is argued that a diffuse electron population that is uniformly distributed over the entire 
supernova remnant of Cassiopia A can contribute to radio and infrared fiuxes. Also in ref. [38, 39] it is shown that 
in the infrared region of Cen A the best fitted photon spectral index lies around 0.89 ± 0.25 and 0.84 ± 0.18 which 
further strengthen the non-thermal synchrotron origin of the infrared emission. 

It has been suggested that both radio and X-ray emission in the inner jet of Cen A are produced by relativistic 
electrons with Lorentz factor 7 ~ 8 x 10^ in 60 iiG magnetic field 40]. With our fitted parameters we found that the 
magnetic field lies very close to 60 nG as shown in Eg.© and at the same time the electron Lorentz factor at the 
break energy is 7 ~ 2.3 x 10^. The maximum observed photon energy for the above Cen A parameters is given by 

<™a. = ^^^^ ^ 1-65 GeV, (6) 

where 77 ^ 1. If the second peak of nuclear SED is due to electron synchrotron emission, the maximum observed 
photon energy from this process is 1.65 GeV and the synchrotron photon spectrum is given by 

(7) 

{ \EoJ ^7 - ^l,b 

where 

and Ey^b — 150 keV is the break in photon energy corresponding to the break in electron energy. Inclusion of the 
first break at far-infrared energy in Eq.0, the condition E^ < E^^i, will be replaced by 4 x 10"'' eV" < E~^ < E^^b- 
During the period 1991-1995, the observation of Cen A by Compton Gamma-Ray Observatory (CGRO) revealed the 
break in the spectral energy distribution vF^ around 150 keV with a maximum flux of about ^ 10~^ erg/cm^ /s[36|. 
By considering a mean power law index ~ 1.85 the observed luminosity in the energy range 40 keV to 1200 keV is 
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E^ < Ey^b 
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\CeEoJ 


Ej > E~f^b 



j^obs _ 3 jg ^ W^^erg/slS&l where the differential photon spectrum can be given by 



7 

-A 



with A^ ~ 5.87 x 1Q~^ / keV/ cm'^ / s . But it is observed that below 150 keV the power index A ^ 1.8 — 2.0 and for 
Ej > E^^b it is 2.4 ± 0.28 . The observed spectrum is related to the synchrotron spectrum at the source as 



obs ,-r^4/ A ,obs\2 

dE^ J dl \ dEy 



dN^Y"" 4r4(At°^^)2 _^^^^ fdN^ 



e--- [^] , (10) 

sy7i 



where dz = 3.8 Mpc is the distance of Cen A from us, and Tint is the optical depth of photon in the jet. The term 
g-Tint [^ ^j^g above equation corresponds to the decrease in observed spectrum due to 77 interaction in the photon 
background of the jet. In this energy range we have observed that Tmt < 10^^ and does not reduce the observed 
photon spectrum'41] which is shown in Figure 1 for different observed luminosities. So for low energy, one can neglect 
the contributions due to 77 interactions. Here we use a — 2.7 as seen in the diffuse UHECR flux just before the 
Greisen-Zatsepin-Kuzmin cut off assuming that like proton, electron also follows the same power spectrum. So in 
Eq.(l7]), a = 2.7 is consistent with the observed with the observed photon spectrum as described in Eq.®. Comparison 
of Eq.dlD with Eq.^ gives A^ ~ 5.4 x 10^^ /keV/cm'^/s. 

In the Inverse Compton (IC) scattering process, the energy of an IC photon is very high and for 40 keV nfinimum 
observed photon energy the electron Lorentz factor is 7 '^ 6.2 x 10^. So for IC to be operative in this region, the photon 
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FIG. 1: We have shown the internal photon optical depth r^-y as a function of photon energy E^ in the nuclear region of 
the jet for three different luminosities. These are (a) continuous line is for L°''° = 3.78 x 10^^ erg/ s, (b) Dashed line is for 
1°^" = 1.18 X 10''^ erg/s and (c) Dot-dashed line is for 1°^" = 2.0 x 10^^ erg/s. 



energy has to be really very small and for this reason we do not consider the contribution from IC process |32ll42l]. The 
combined OSSE (0.05 - 4 MeV) ^S], COMPTElHI and EGRET (0.1 - IGeV) jQ) data is fitted by doubly broken 
power law with two break energies (~ 150 keV and 16.7l^g'3 MeV). The spectral indices are A ~ 1.74 for E^ < 150 
keV, A ~ 2.3 for 150fceF < E^ < WJtlH MeV, otherwise A ~ 3.3 for E^ > U-TtlH MeV. But the second break at 
Ej > 16.7l^g| has large uncertainty and also the spectrum is almost smooth around and beyond this energy. So in 
the present analysis we do not take into account for this break and assume that the power index is 2.35 well fitted. 

We have shown in Eq.Q that the electron synchrotron process contributes up to maximum observed photon energy 
~ 1.65 GeV. Also in the nuclear region r ~ lO.Spc, the inelastic collision of Fermi accelerated protons with the 
background protons will produce 7-rays due to decay of neutral pion sllSi. I n this case the minimum observed energy 
of proton is E"''" ~ 5.2 GeF of which about 18% goes to a photon|44 l45l| which corresponds to E"''" ~ 0.93 GeF. 
So for £'"'"' > 0.93 GeV, pp collision will also produce 7-rays and the projectile high energy protons (due to Fermi 
acceleration) do not have a break energy (the only break in proton energy will come from the proton synchrotron 
emission which is given by Ep^i, ~ 2.8 x 10^^ eV). In the nuclear region the same Fermi accelerated protons colliding 
with the background photon should also produce photons through delta resonance and its flux can also compete 
with the one from pp process. By considering the highest detected energy bin found by the Energetic Gamma-Ray 
Experiment telescope (EGRET) at a break energy ~ 200 MeV [g, [ij|, we observed that the p^ process starts to 
contribute at energy E^^rnin — 1-4 GeV. But the flux ratio oi pp to p"f is about ~ 3.13 x lO^'' at 1.4 GeV photon 
energy which increases by a very small amount for higher energies which shows that p^ process is negligible and does 
contribute to the very high energy gamma rays. 

In the observed photon energy range 0.93 GeV to 1.65 GeV both electron synchrotron and pp interaction will 
contribute and the pp contribution can be given by 
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where 
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5 = 0.18 and Up is the proton density in the jet and CTpp ~ 4.5 x 10 



-26cm2. 



(12) 
The energy normahzation we have used 



here is £^0 = 1 GeV. Around r ~ 10.8 pc from the central black hole the proton density can be in the range of 10 /cm 
to 10^/cTO'^ ref . [ij, [S^] . Here we consider an intermediate value of Up '^ 10^/cm'^ and the optical depth of the process 
is Tpp = rupapp ~ 0.15, which shows that about 15% of the protons will interact to produce photons. In this energy 
range the photon power spectrum from both electron synchrotron and pp are proportional to E^ and E~°' 

respectively. Also it is important to note that we have not observed any break in photon energy around 0.93 GeV and 
beyond and on top of that the photon power spectrum is well fitted with the index 2.4 ± 0.28 in the energy range 30 
MeV to 10 GeV [g. The EGRET has measured the average flux above 100 MeV (100 MeV to 10 GeV) from Gen A 
which is N-y = (13.6±2.5) x \Q~^ photons / cm? / s and for A = 2.35 the luminosity obtained is L"'"' = 1.19 x 10^^ erg/s. 
As we have discussed earlier, from 100 MeV to 0.93 GeV, the luminosity is only due to the electron synchrotron 
radiation which is given as L^ ~ 9.24 x lO'*" erg/s. On the other hand in the energy range 0.93 GeV to 1.65 GeV, 
the luminosity from electron synchrotron is L^ ~ 1.85 x IQ'^^ erg/s and from pp is L^ ~ 3.23 x 10'^^ erg/s. Finally 
the luminosity in the range 1.65 GeV to 10 GeV from pp process is L^ 



the parameters Up and Ap so that photon spectrum from pp interaction will overlap with the one from electron 



synchrotron. But once we fix the value of the parameter Up in the nuclear region (we take Up c^ 
has to be adjusted. The observed photon spectrum in the energy range 0.93 GeV to 1.65 GeV is 



4.69 X 10'^^ erg/s. So we have to adjust 
L the O] 

10^/cm^), only Ap 
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The fitted value of Ap is given by Ap = 1.8 x XQ^ /GeV/cm? / s. Now we shall turn to the high energy range i.e. 
above 250 GeV in which many collaborations have measured the photon flux and also the spectral index. As we have 
discussed p^ process is negligible compared to the pp process, so only pp process will contribute to the very high 
energy gamma-ray production in Cen A. For E^ > 250 GeV only pp process will contribute. But before that we would 
like to estimate the number of events from PAG which is shown in the next section. 

III. COSMIC RAY EVENTS FROM PIERRE AUGER OBSERVATORY 

The Pierre Auger Observatory (PAO ) ha d reported a correlation between ultra high energy cosmic rays (UHECR) 
and nearby AGN within ~75 Mpc[4y, |43|- Roughly 10 cosmic ray events are concentrated around the Centaurus 
direction, a region with a high density of AGNs and two of these events fall within 3 degrees from Cen A suggesting 
that it could be the first evidence of the UHECR source. Although the composition of these CR events are unknown, 
we assume them to be protons |16ll2l|. 

In general the intergalactic magnetic field plays a major role in deflecting the CR protons from their original 
direction. But the deflection angle AO is inversely proprotional to the CR proton energy and for a 57 EeV proton 
it is less than a degree, if we take the intergalactic magnetic field to be of order 10~^ G and for higher energy CR 
protons, the deflection will be still smaller. So the propagation of CR protons from Cen A with energy > 57 EeV are 
almost ballistic and one can neglect the effect of magnetic field on it. 

We estimate the expected number of events from Auger array. For a point source, the integrated exposure of PAO 
is 5 = 9000/7rfcTO^ yr and the integral flux of UHECR from PAO can be expressed as 



N(> Ermn) 



^ Events 
Exposure ' 



(14) 



where E„ii,i = 57 EeV. One has to also consider the relative exposure uj{6) for angle of declination 6. For Cen A, 
S = 47° and the corresponding value of w((5) ~ 0.64[14]. The time duration for the collection of data by PAO is about 
15/4 yr between 1st January 2004 and August 2007. So from the above equation we obtain 



,/ rn 9 X lQ^km?yr ^ . ^ , ,„ 15 _ 

# Events = A(> A„ji„)a;(d) -—e 
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(15) 



where e "^p^ is the survival probability of proton due to interaction with the protons in the jet background. Also the 
integral flux is related to the proton normalization constant Ap as 
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(16) 



Putting Ea. (|T6|) in Ea. (fT5|) and simplifying we obtain 

# Events = 5.75 x 10^ x Apo, (17) 

where we have defined Ap = Apo GeV^^ cm^^ s^^ . We can obtain the value of Apo from different observa- 
tions/Experiments which are discussed in the next section and using these values we can calculate the expected 
number of events in these observations. 

IV. OBSERVATION OF CEN A ABOVE 250 GEV 

As discussed in the introduction, Cen A has been observed in almost all wavelengths. Following are the observations 
of Cen A in Very High Energy 7-rays (> 250 GeV) |48-50] . From each observations we have used their given integral 
flux and calculate the photon normalization constant A^ as given in Eq. ^ and after that comparing with the observed 
spectrum oi pp process calculate Ap. 

1. The H.E.S.S. Collaboration operates an array of four large imaging Cherenkov telescopes to detect Very High 
Energy (VHE) 7-rays in located in Southern hemisphere in Namibia. Cen A was observed between April 2004 
and July 2008 with total live time of 115.0 hours. The measured integral flux above 250 GeV energy threshold 
was (1.56 ± 0.67) x 10~^^ cm~^ s~^ which was around the central region of the galaxy. Also the measured 
differential photon spectrum is well described by a power-law with the index ~ 2.73 and the normalization con- 
stant A^ ~ 2.45 X 10^1^ /TeF/cmVsiT]. This gives the luminosity L^ - 2.58 x lO^^ erg/s and Apo =0.18x10-3. 

2. In March 1997, with the Durham Mk6 telescope, Cen A was observed for 6.75 hours with a 3a flux upper limit 
above 300 GeV of 5.2 x 10"" cm-^ s-^^]. This gives A^ ~ 1.14 x 10^" /TeV/cm^/s. 

3. The JANZOS group observed Cen A from New Zealand during April 1988 and June 1989 for 56.9 hours and 
reported a 95% confidence level upper limit on the flux above 1 TeV of 3.74 x 10^^^ cm'"^ s~^ [10] which gives 
A^ ~ 1.14 X 10-" /TeV/cm^/s. 

4. In March and April 1999, Cen A was observed by the CONGAROO 3.8 m telescope for 45 hours with a 3a flux 
upper limit above 1.5 TeV of 5.5 x 10^^^ cm-^s-^ for a point source at the core of the galaxy [11[ and this 
gave A^ ~ 1.86 x 10"" /TeV/cm'^/s. 

5. In March and April of 2004 further observations were made to Cen A with three 10 m telescopes of the 
CONGAROO-IH array for 10.6 hours. For the core region of Cen A with a 2a flux upper limit above 424 GeV 
was 4.9 x 10-12 cm-^ s'^ ^ and we obtain A^ ~ 1.94 x 10-^^ /TeV/cm^/s. 

In all of the above VHE 7-ray observations we have fitted the differential photon spectrum with the power index 
a = 2.7[3, ImI l52| and calculated their normalization constant A^. Comparison of A^ with the observed photon 
spectrum from the respective process gives the values of Ap which are given in table 1. Using this Ap we calculate 
the number of events from these different observations. 

It is observed that the H.E.S.S. collaboration results give less than one event during the observation period of 
15/4 years. On the other hand during the same observation period the Durham Mk6 telescope, CANGAROO-IH, 
JANZOS and CANGAROO observations could have given respectively 5.06, 0.86, 16.61 and 8.26 events in PAO due 
to pp interaction in Cen A. So our analysis shows that many UHECR protons from Cen A can be observed by PAO. 

V. CONCLUSIONS 

The Cen A is observed in almost all the energy band and its nuclear SED shows two peaks of which one is around 
150 keV and another is in the far-infrared band. In this work we assumed that the peak around 150 keV is due to the 
electron synchrotron radiation and adjust the parameters F, ^^ , Ai"*"* and fes so that the observed photon energy 
E"^" is 150 keV. In this fitting we have F = 4.2 and /Sl"^" = 1 yr correspond to the fact that the jet is in the nuclear 
region r ^ 10.8 pc and it is mildly relativistic which is being observed in Cen A. The magnetic field obtained around 



TABLE I: Observation of Cen A by different Experiments with their respective lower observed photon energy limits are shown 
here. We also show their corresponding proton normalization constant and the expected number of CR proton events. 



Experiment [rof.] 


E-, 


Apo 


# Events 


H.E.S.S.f7] 


> 250 GeV 


0.18 X 10"^ 


0.1 


DurhamfS] 


> 300 GeV 


8.81 X 10~^ 


5.06 


CANGAROO-III [9] 


> 424 GeV 


1.5 X 10"^ 


0.86 


JANZOS [10] 


> 1 TeV 


28.9 X 10"^ 


16.61 


CANGAROO fll] 


> l.STeV 


14.38 X 10"^ 


8.26 



the nuclear region is about 180 ij,G[AO]. We found that the IC process does not contribute to the observed photons 
in the GeV- TeV range. In our model the maximum energy up to which the electron synchrotron contributes is 1.65 
GeV and the other process which contributes to high energy (1 to 10 GeV) to very high energy (> 250 GeV) is the 
pp process. Unfortunately the P7 process is very small compared to the pp process. 

We compared the observed photon spectrum from different observations by taking the spectral index a = 2.7 with 
the calculated one from the pp process, which determines the normalization constant Ap for the cosmic ray (CR) 
protons. We use the Ap from these observations to calculate the number of events from PAO and found that indeed 
many CR protons above 57 EeV can arrive on Earth from Cen A. In a recent paper in the context of Cassiopia A, a 
young supernova remnant, Araya and Cui[37] have also independently come to the same conclusion that synchrotron 
emission from the electron can account for the radio to X-ray emissions and additional hadronic component is needed 
to explain the GeV emission thus providing evidence for the production of cosmic rays. 

We are grateful to S. Nagataki and Y. Y. Keum for many useful discussions. This work is partially supported by 
DGAPA-UNAM (Mexico) Project No. IN101409 and Conacyt Project No. 103520 (SS) and by NAS A NNX09AO94G 
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